DNA translocation through nanopores with salt gradients: The role of osmotic flow 
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Recent experiments of translocation of double stranded DNA through nanopores [M. Wanunu et al. Nature 
Nanotech. 5, 160 (2010)] reveal that the DNA capture rate can be significantly influenced by a salt gradient 
across the pore. We show that osmotic flow combined with electrophoretic effects can quantitatively explain the 
experimental data on the salt-gradient dependence of the capture rate. 



Translocation through solid-state nanopores holds the po- 
tential to be a fast commercial method for macromolecular 
characterization and sequencing, such as for long, unlabelled 
single-stranded and double-stranded DNA molecules fl-fstl- 
Clearly, high throughput and time resolution — effected by 
enhanced capture rate as well as translocation times respec- 
tively — is a necessary precondition for the process to be 
commercially viable. Although the capture rate or transloca- 
tion time can be increased by manipulation of the temperature, 
salt concentration, electric field strength and viscosity |4|. the 
increase of one is usually accompanied by a decrease of the 
other [4]. Recently however, Wanunu et al. jsl showed that 
it is possible to increase the capture rate without decreasing 
the translocation time by using a forward salt concentration 
gradient across the pore. The large increase in capture rate 
as a function of the imposed salt concentration gradient was 
qualitatively explained by the increase in the electrophoretic 
motion of DNA towards the pore as a function of salt asym- 
metry: A constant current of ions is flowing through the pore, 
creating a long range electric field which acts as a funnel for 
the ions and the polymers towards the pore js]- 



These studies U-iZl], and some others on similar sys- 
tems lil- lioll . have focused on electroosmotic and/or elec- 
trophoretic phenomena. Electroosmotic phenomena describe 
flow of liquids with a net mobile charge under applied electric 
fields, while electrophoretic effects relates to the movement 
of charged polymers in an electric field. [These terminologies 
are further discussed in Sec. I of the Supplementary Online 
Material]. It is worthwhile to note at this point that electroos- 
motic flow in the experiments of Wanunu et al. was found 
to be in the opposite direction of the observed DNA translo- 
cation (see the supplementary material of Ref. jsfl), imply- 
ing that electroosmotic flow cannot explain the observed en- 
hanced capture rate. In fact, in the presence of an imposed salt 
concentration gradient across the pore, there is an additional 
mechanism at play, namely the capillary osmosis process ifTll - 
1311 . i.e., the flow of water from a lower osmotic pressure 
(cis) side to a higher osmotic pressure (trans) side through 
the pore. The effects of this osmotic flow on the DNA cap- 
ture rate has been missing in the theoretical analysis so far. 
In this Letter we show that the osmotic flow is a key ingredi- 
ent to understand the experiments of Wanunu et al. jstl; their 
results can be quantitatively explained with osmotic flow and 
electrophoretic effects 11211 . Note also that the full range of 



dynamical mechanisms affecting DNA capture in the experi- 
mental setup of Wanunu et al. are discussed in Sec. I of the 
Supplementary Online Material. 

The osmotic flow of water is driven by a pressure gradient 
anti-parallel to the salt concentration gradient. The reservoirs 
are kept at a constant pressure, such that a chemical poten- 
tial gradient is present accross the pore for both the ions and 
water, causing flow of ions down the salt concentration gradi- 
ent and water up the salt concentration gradient. However, a 
pressure gradient inside the pore and a correspondig net flow 
of the liquid (ions plus water) will only be present if the ions 
are net depleted from (or net attracted to) the pore 11411 . In 
the textbook example of osmosis through a semi-permeable 
membrane ions are completely restricted from entering the 
pore due to steric repulsion 11411 . However, also when the 
restriction is only partial, e.g. due to wall-ion interactions in 
the nm vicitiity of the pore walls, an osmotic flow develops 
fill [12I. [Isl Ea). Water and ions confined in a nano pore can 
behave very differently from the same bulk system 117142311 . 
Both water and ions will be influenced by the pore walls, lead- 
ing to attraction or depletion of ions and/or water. To capture 
such a behavior with a simple continuum model we introduce 
one length scale describing the interaction of the ions with the 
pore wall (the depletion length), and one length scale quanti- 
fying the slip of water flow at the pore wall (the slip length). 
To quantitatively describe the experimental data of Wanunu et 
al. lHll, the ion-wall interactions is found to be repulsive, in 



agreement with simulations Il6l 12211 and theories 12111241] of 
kosmotropic (hydrophilic) ions near low dielectric surfaces. 
Also, the flow of water is found to be in the slip-flow regime, 
in agreement with experimental studies of flow at smooth sur- 
faces I25II . In the experimental system studied by Wanunu 



et al. U we find the osmotic flow to provide the dominant 
contribution to the enhanced capture rate for weak salt gradi- 
ents. In the same experiment [SQ electroosmosis was found 
be a weak effect, reflected in a two/three orders of magnitude 
lower capture rate of neutral PEG compared to charged DNA. 
Also the measured purely Nernstian behavior of the diffusion 
potential suggests nearly neutral pore walls. 

The geometry we study, similar to the experimental setup 
of Ref. fsll, is shown in Fig. [T] Two reservoirs at constant 
pressure Pq with salt concentration Ct (trans side) and Cc 
(cis side) are separated by an impermeable solid membrane 
of thickness L. A cylindrical pore of diameter d connects the 
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two reservoirs. The two electrolytes are composed of mono- 
valent ions of concentrations Ca, and J2a=± ~ ^- The 
solvent (water) is modeled as a continuum with dielectric con- 
stant e = 80, and viscosity i] at temperature T. The Debye 
screening lengths n^J^, are defined as K^y^ = AnCc/tl^e^ / e, 
where /3^^ = ksT, ks is the Boltzmann constant and e is 
the elementary charge. Due to the preference of ions to be 
solvated in bulk water, they feel a repulsive potential U (p) 
from the pore walls 1119112 ill , where p is the radial coordinate 
around the cylindrical axis inside the pore, measured from the 
center of the pore. We model such interactions with a region 
£ next to the pore walls depleted of ions (see Fig. [TJ. The 
part of the pore accessible to ions is described by the diame- 
ter a = [d — 2£). The polymers (DNA) are located in the cis 
chamber, and the electric field is applied from the trans to cis 
side, driving DNA (with negative charge) from the cis to the 
trans reservoir 

The pore is assumed to be neutral, i.e. J2a QaCa{p, z) = 
istl, where q± = ±e, which is a good approximation for nearly 
neutral pore walls and low induced charge within the pore. 

For the system studied here, the Reynolds number is very 
small, and the flow can be described by the steady Stokes 
equation combined with incompressibility of the liquid: 

7?V2v(P,z) = VP(p,z) + ^c,(p,z)VC/(p); (1) 

a 

V-v = 0. (2) 

In steady state the dynamics of the ions are described by the 
time-independent Nernst- Planck equations: 



+V-(c„(p,z)/3[VC/(p)-<z„E(z)])) (3) 
+V- [c,(p,z)v(p,z)] =0. 

This is equivalent to conservation of particle current, with J a 
the current density. Da the diffusion coefficient of ion type a 
and E(z) the local electric field. Assuming fast equilibration 
of the concentration and the pressure in the radial direction 
(/5 • Jq = and p • v = 0) we get from Eqs. ([B and (O HIl] 
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Ca{p,z) R 

r,y^Vp{p) = 0, 
which from Eq. (|TJ gives 

^^^(^'")^{^BmCo(.) pta/l. 

If we further assume that the ion density changes linearly ac- 
cross the pore (which follows from Eq. (O when diffusion 
dominates over convection) we get 
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FIG. 1: Schematic of the pore geometry showing a membrane of 
thickness L connecting two salt reservoirs with salt concentrations 
Cc(cis) and Ct (trans) by a pore of diameter d. In the cis reservoir 
there is a bulk DNA concentration of Cdna, and a voltage difference 
V is applied accross the system. The salt ions are depleted within a 
layer £ from the pore walls. There is a liquid velocity profile v^{p) 
in the z direction, which varies with the radial coordinate p. There is 
slip of the flow at the pore walls described by the slip length b. 



Since the pore is in the nanometer regime a continuum treat- 
ment of the fluid dynamics may not be accurate. The Knudsen 
number for the system is Kn = S/d ~ 0.1, where (5 « 3 A is 
the intermolecular spacing for water ll26ll , indicating that we 
are in the slip-flow regime (0.01 < Kn < 0.1), such that 



v,id/2) 



dv^ip) 



dp 



(8) 



p=d/2 



where b is the slip length (see FiglT]!. The flow can now be ob- 
tained by integrating Eq. ^ twice making use of Eq. ([8]l. The 
resulting area-averaged velocity of the flow inside the pore is 



Vo 



kBT{Ct~Cc)aodHl + 8b/d) 
L 3277 



(9) 



where we have introduced the osmotic reflection coefficient 

""^^^ i+itb/d) ^^^^Z'^) + ^ ~ ■ ^^^^ 

For a = (ctq = 1), i.e. ions are totally depleted from the 
pore, we recover the standard slip-modified Poiseuille flow 
due to osmosis through a semipermeable membrane 12711 . If 
we set the slip length to zero, we recover the result of An- 
derson and Malone for leaky membranes llT4ll . however with 
an effective solute radius £. With £ = (no ionic depletion) 
Eq. ( [Tol l gives ctq = (no flow), showing that ion depletion 
is crucial for the present analysis. From Eqs. ([T]i and (|2]i the 
osmotic flow at a radial distance r ^ d from the pore can now 
be approximated as 



.Vod"^ 



(11) 
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where f is the radial unit vector, pointing outward from the 
pore mouth. 

In a steady state and using conservation of charge current 
(Eq. (|3]l), the electric field on the cis side (for \r\ ^ d) can be 
approximated as fsf] : 
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E(r) = 



^ Cpa^V 



(12) 



where Cp — (Ct + Cc) /2 is the ion concentration inside the 
accessible part of the pore, and V/L is the strength of the 
applied i?-field in the pore. The drift of charged polymers in 
an electric field is described by electrophoresis 128 1 



E(r), 



(13) 



where 0dna is the surface potential of DNA, and n is the 
electrophoretic mobihty. 

To get an estimate of the number of polymers that translo- 
cate through the pore per second, we calculate the flux of DNA 
generated by the combination of electrophoretic effects and 
osmosis. By conservation of DNA particle current, we get the 
capture rate per bulk DNA concentration 



Rr 



27rr^ VEp(r) + vos('') 



(14) 



independent of r. Flow towards the pore is antiparallel to f 
(see Fig. [T]i, and therefore a negative sign in Eq. ( fT4] i ap- 
pears such that Rc > for translocation from the cis to trans 
reservoir Combining Eqs. ( fTTT i. ( fT3] l and (fl4] i we find: 



Rcix) = Rc{l) 
where x ~ Ct/Cc and 
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(16) 
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Note that the result does not depend on DNA length [5l|, and 
that fc = (or a; = 1) describes electrophoretic effects alone. 

In Fig. 12] we plot the predictions of Eq. ( fTSl l as a func- 
tion of salt asymmetry x for different values of the dimen- 
sionless parameter k (Eq. ([TT])). As the value of k increases 
the predictions start to deviate from the capture rate due to 
electrophoretic effects alone (straight line, k = 0). The flow 
due to osmosis varies inversely with x (i.e. linear in Cc), and 
will therefore saturate for large x, while the flow due to elec- 
trophoretic effects has a linear dependence on x with slope 
1/2. With k in Eq. ( fTSl l as a free parameter, we find an excel- 
lent fit to the experimental data of Ref. |Q] for fc « 7.5 and 
X < 5. Having noted that the calculations presented in this 
work are valid to first order in the salt concentration gradient 
(Ct — Cc)/L, in Fig. |2]we focus on the regime where x < 5, 
as these data points are all for a 1 M salt solution in the trans 
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FIG. 2: Capture rate (Eq. (TS}) as a function of salt asymmetry, for 
different values of the dimensionless parameter k (see Eq. l|15|l). The 
points are experimental mesurements of Ref. (J. 



chamber, and represent about 75% of the available data. Note 
in this context that the experimental data of Ref. jstl also con- 
tain four additional data points with larger salt concentrations 
in the trans chamber (2 M and 4 M), and larger values of x. 
These data points lie outside the region where we expect our 
model to be valid. 

To further compare our predictions with the recent experi- 
mental measurements of DNA translocation in salt gradients 
ijstl, we fix the system parameters to the experimental values 
(for DNA length N = 400 bp and N = 2000 bp); Ct = 1 M, 
d 3.5 nm, V = 300 mV and L = 20 nm. For the DNA 
electrophoretic mobility we use ji = — lO^^m^s^^V^^ ll29ll 
which for water with viscosity rj ~ 1 mPa ■ s (at T = 20°C) 
gives from Eq. (fTsl l /3e(/)DNA ~ —0.55. The only free pa- 
rameters are the depletion length i and the slip length h. For 
depletion lengths £ = 0.3 to 0.6 nm, which one expects due to 
finite-size effects of (hydrated) ions like K+ and Cl~ l24ll and 
image charge effects 12111 . one finds from k ~ 7.5 that & = 4 
to 7 nm (see SOM for details), in reasonable agreement with 
measured values of the slip length at smooth surfaces ll25ll . In 
the Supplementary Online Material we plot combinations of 
I and h corresponding to different values of k, as well as the 
osmotic flow profile Vz{p) (see Sec. II). 

The ions are assumed to be depleted from the pore walls, 
an assumtion which is based on numerous theoretical Ell |23[ 
2411 . simulation 1122 , and experimental studies Il34l35ll . 

that find nonpolarizable ions to be repelled by an interface 
between water and a low dielectric material such as Silicon 
Nitride (e = 7). This behavior can be modified due to surface 
chemistry, such as dangling atoms Bill , surface charge and 
affinity for water [22J. Dangling atoms lead to binding of ions 
to the surface, which can result in current rectification 131 1. 
Unless these effects are very strong, the ions are generally all 
over depleted from neutral low dielectric interfaces. 

To conclude, having approximated the ion-wall interactions 
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due to image charges and water structure by an effective de- 
pletion length we show that the experimental data of Wa- 
nunu et al. for DNA translocation in salt gradients can be 
explained by a combination of electrophoretic effects and os- 
mosis. To account for the different behavior of water on the 
nanoscale we have also introduced hydrodynamic slip at the 
pore walls, which enhances the flow due to osmosis. With 
resonable values for both the slip length and the ion depletion 
length, we find quantitative agreement between theory and ex- 
perimental measurements. 

Throughout the calculation we have focused on the diffu- 
sion limited regime, and do not take into account the free- 
energy barrier felt by the polymers when entering the pore, 
yet we are able to quantitatively reproduce the relative cap- 
ture rate enhancement data in the barrier-limited regime. This 
is most likely a signature of the fact that the barrier height is 
nearly constant as a function of salt asymmetry. We do expect 
that the main physics reported here, namely the role of the os- 
motic flow, explains the enhanced relative capture rate in the 
barrier limited regime; however, including the barrier in our 
analysis remains a significant challenge. 

We have also assumed the ion density inside the accessi- 
ble part of the pore to be equal to the average of the salt 
concentration in the two reservoirs. This assumption is sup- 
ported by the current-voltage relations measured by Wanunu 
et al. for different salt concentrations in the cis chamber (for 
Ct = IM, Cc = 0.2Mto IM), see supplementary information 
of Ref. IS\. The calculations presented in this work are valid 
to first order in the salt concentration gradient (Ct — Cc)/L, 
and the results of the measurements by Wanunu et al. with 
higher salt concentration in the trans reservoir (x > 5), is out- 
side the region where this model is expected to be valid. 

Putting things in perspective, translocation of DNA through 
nanopores is a complicated problem due to its many aspects, 
ranging from properties of water in confinement to compli- 
cated structures of the translocating molecules and their inter- 
actions. However, to understand the recently found increase 
in capture rate as a function of salt concentration asymme- 
try, it seems that a detailed description of DNA molecules is 
not needed, since the main mechanism is the enhanced attrac- 
tion of DNA molecules towards the pore. This attraction is 
here shown to be made up of two main contributions: elec- 
trophoretic effects and osmotic flow. The capture rate with salt 
gradients due to electrophoretic effects ISl |7[] and combined 
with electro-osmosis has been described before, however 
the role of (diffusio) osmosis has not been previously dis- 
cussed. To understand the osmotic flow it is crucial to account 
for the repulsive interaction between kosmotropic ions and a 
neutral nonpolar wall. We expect that the main physics is cap- 
tured by introducing a layer near the pore wall depleted of 
ions. This also means that the osmotic flow is ion specific, and 
will be reversed when using salt particles that are net attracted 
to the pore wall. Finally, our analysis shows that osmosis can- 
not be ignored for nanopores in the presence of salt gradients, 
even though salt is able to flow through the pore. 
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SUPPLEMENTARY ONLINE MATERIAL 



In the main text we discuss the dominant contributions to 



enhanced DNA translocation induced by a salt gradient in the 
recent experiment by Wanunu et al. yj]. In this Supplemen- 
tary Online Material we also discuss other relevant dynamical 
processes that can drive DNA towards a nanopore in the pres- 
ence of both an external electric field and a salt concentration 
gradient accross the pore. We also show plots of the induced 
osmotic flow inside the pore, and relevant combinations of the 
ion depletion length £ and the slip length h (see main text for 
details). 



I. RELEVANT DYNAMICAL PHENOMENA FOR DNA 
TRANSLOCATION IN THE PRESENCE OF SALT 
GRADIENT AND A DRIVING VOLTAGE 

For a narrow pore in the presence of an external electric 
field, there are two dynamical phenomena that can drive DNA 
towards/away from the pore: 

(i) electroosmosis, and 

(ii) electrophoresis. 

These two mechanisms can be put under the broader group 
of electrokinetic effects, which describe the motion of fluids 
and charged solutes under forces of electrostatic origin. Phe- 
nomenon (i) is the movement of a fluid medium containing a 
net mobile electric charge driven by an external electric force; 
for global neutrality reasons this requires the presence of fixed 
charged walls, as the screening cloud in the fluid constitutes 
the net mobile charge that couples the electric field to the fluid 
motion. Phenomenon (ii) concerns the motion of a (often, 
confined) charged polymer, such as DNA or RNA, under an 
applied electric field, e.g., gel electrophoresis, used for seg- 
regating DNA by their lengths; it describes the motion of a 
DNA molecule through the pores of a gel under an external 
electric field. 

In the presence of a salt gradient two further dynamical 
mechanisms play a role for the DNA capture by the pore: 

(iii) osmosis, and 

(iv) diffusiophoresis. 

Osmosis (iii) is the flow of water from a low osmotic pres- 
sure to the high osmotic pressure side, i.e., from the cis 
to the trans in the experiments of Wanunu et al. Osmosis 
through membranes that are permeable to ions is often called 
capillary-osmosis or diffusio-osmosis |f3l, to distinguish the 
mechanism from classical osmosis through a semipermeable 
membrane. Similarly, associated with the concentration gra- 
dient of solutes is diffusiophoresis (iv), which describes the 
movement of macromolecules in a solute concentration gra- 
dient fst]. For both diffusiophoresis and (diffusio) osmosis, 
the ions (solutes) must either be attracted to or repelled by 
the wall/macromolecule. This repulsion/attraction leads to a 



lateral pressure gradient in the region where the ions interact 
with the wall/macromolecule, driving fluid flow or movement 
of the macromolecule. In the main text we discuss osmosis 
(iii) and electrophoresis (ii), which we find to be the dominant 
contributions to the enhanced DNA capture rate for the partic- 
ular system studied by Wanunu et al. ^ i.e., electroosmosis 
is subdominant. Below we argue that diffusiophoresis (iv) can 
also be neglected in relation to electrophoresis. 

The flow due to electroosmosis as driven by an electric field 
Ep in the pore has the exact same form as the flow due to 
electrophoresis (with a change of sign) 101 



VEO 



Ami 



(18) 



so when both the polymer and the pore have surface charge 
(surface potential) with the same sign, electrophoresis and 
electroosmosis pull the polymer in opposite directions. In 
the experiment by Wanunu et al. lU the electroosmotic flow 
was found to be negligible compared to electrophoresis, as 
the pore walls was found to be nearly neutral. The diffu- 
siophoretic flow, which relates to diffusioosmosis as elec- 
trophoresis relates to electroosmosis, is induced by the salt 
gradient in the cis reservoir, and pushes DNA towards the pore 
(as ions are net attracted to DNA molecules). By conservation 
of particle current 



J = -i:>VC(r), 



(19) 



the gradient of the salt concentration in the cis reservoir is 



VCc(r) = 
and the corresponding flow is 



.{Ct-Cc)d^ 



vdp(»') = £'DpVlnCc(r) = -D^pf 



(^-1) 
8ir2 ' 



(20) 



(21) 



where Ddp is the diffusio-phoretic mobility. The diffusio- 
phoretic flow has the same dependence on x as electrophore- 
sis, and will therefore add to the term linear in x in the relative 
capture rate (see Eq. (15) in the main text). The ratio between 
the electrophoretic flow and the diffusiophoretic flow is (from 
Eq. iTH and Eqs. (12) and (13) in the main text) 



|vDp(r)| _ 2Dj^p{x-l)d^ 



IvepWI 



nVa'^x 



0.1 



(22) 



with the electrophoretic mobility of DNA i^i = 
-lO-^m^s-iy-i, V = 300 mV and Ddp « 2 • IQ-^^ 
m^/s fst]. This is about one order of magnitude lower than 
electrophoresis and is therefore omitted in the analysis in the 
main text. 



II. OSMOTIC FLOW PROFILE 

In Fig. Oa) we plot the different combinations of £ and 
b for different values of k, showing resonable combinations 
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FIG. 3: (a) Combinations of the slip length b and ion wall depletion length £ for different values of k (see Eq. (17) in the main text), with all 
other values fixed to match the experimental values of Ref. 1 5]. (b) The liquid flow profile inside the pore as a function of the radial coordinate 
p for different combinations of the slip length and depletion length with k = 7.5 and x = 5, all other parameters are fixed to match the 
experiments of Ref. ^ (see main text for details). 




C [M] 

FIG. 4: Capture rate (Eq. (15) in main text) as a function of salt 
concentration in the cis reservoir, with Ct = 1 M, for different values 
of the dimensionless parameter k (see Eq. (15) in the main text). The 
points are experimental mesurements of Ref. |@| 



of i and b for the fc-regime of interest (see main text for de- 
tails). The slip length of surfaces are typically in the range 
— 30 nm depending among others on the contact angle and 
the smoothness of the surface |i5i|. Similarly, in Fig. [3jb) the 
liquid flow profile inside the pore is plotted for different com- 
binations of the depletion length £ and the slip length b, with 
k = 7.5, X = 5, and all other parameters fixed to match the 
experiment of Ref. |^ . When the slip length is larger than 



the pore diameter, b > d, the flow profile becomes almost uni- 
form, in contrast to the no slip case where the flow is zero near 
the pore walls. When £ is small compared to the pore diam- 
eter a larger slip length is needed to produce a large enough 
flow to explain the experimental measurements. This is be- 
cause a large £ also decreases the electrophoretic motion of 
the DNA molecules, by lowering the ion-acessible area of the 
pore. For the flow profiles in Fig. [3] the Reynolds number is 
Re ~ vd/v < 0.0035, where v = lO^^'m^/s is the kinematic 
viscosity of bulk water. 

In Fig. ID we plot the relative capture rate as a function 
Cc, showing a linear dependence between the capture rate and 
Cc- This is the exact same plot as Figure 2 in the main text, 
however plottet in a different manner This linear dependence 
stems from the osmotic flow as we predict (see Eq. (15) in 
the main text), and not from electrophoresis, suggesting that 
the capture rate has a linear dependence on Cc for weak salt 
gradients, and an inverse dependence for large salt gradients 
(X = Ct/Cc). 
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